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Objectives 

•  Development  of  alternative  approaches  to  the  full  Boltzmann  model  in  transition  regimes. 

•  Applications  to  the  numerical  simulation  of  flight  at  high  altitudes  particularly  the  computa¬ 
tion  of  flows  around  full  airplanes  at  high  altitudes. 


Status  of  effort 

During  the  last  two  years,  the  author  continued  his  efforts  to  reach  the  above  objectives.  In 
particular  he  focused  on  the  following: 

•  Further  numerical  and  theoretical  studies  of  the  coupling  of  Boltzmann  equation  and  its 
hydrodynamic  limits. 

•  Numerical  and  theoretical  analysis  of  the  coupling  of  kinetic  equations  and  their  hydrody¬ 
namic  limits  for  particular  kinetic  models.  (This  gives  more  insight  into  the  coupling  method 
for  the  full  Boltzmann  model) 

•  Numerical  and  mathematical  studies  of  the  coupling  of  BGK  and  ES-BGK  with  their  hydro- 

dynamic  limits.  (The  study  of  BGK  and  ES-BGK  led  the  author  to  replace  in  the  coupling  of 
Boltzmann  equations  and  their  hydrodynamic  limits,  the  Boltzmann  equations  by  ES-BGK 
model.  The  author  plan  now  is  to  show  that  this  resulting  coupling  is  optimal  and  can  be 
used  for  the  numerical  simulation  of  flows  around  full  airplanes  at  high  altitudes).  This  study 
is  subdivided  into  two  parts  -  -  - 

In  the  monoatomic  case,  the  author  is  developing  and  validating  his  coupling  using  Direct 
Simulation  Monte  Carlo  method  for  the  kinetic  model.  He  is  also  planning  on  devel¬ 
oping  a  deterministic  approach  for  the  kinetic  model.  This  will  cut  the  computer  cost 
considerably  and  allows  for  the  computations  around  full  airplanes. 

In  the  polyatomic  case,  the  author  is  developing  and  validating  his  coupling  using  Direct 
Simulation  Monte  Carlo  method  for  the  kinetic  model  when  an  internal  energy  variable 
is  used.  He  will  also  study  a  two  distributions  approach  that  allows  to  bypass  the  use  of 
the  internal  energy  variable  using  Direct  Simulation  Monte  Carlo  method  for  the  kinetic 
model.  Thus  cutting  the  cost  of  computations.  He  will  then  extend  the  deterministic 
method  to  the  two  distributions  approach. 
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1  Introduction 

The  main  problem  of  kinetic  theory  consists  in  deriving  an  evolution  equation  for  the  one-particle 
distribution  function  and  in  solving  mathematical  problems  related  to  this  equation,  so  that  the 
moments  of  the  solution  can  enable  the  calculation  of  the  relevant  macroscopic  variables. 

The  celebrated  Boltzmann  equation  is  certainly  the  most  famous  and  successful  mathematical 
model  to  describe  the  time  and  space  evolution,  at  a  statistical  level,  of  the  position  and  velocity 
distribution  function  of  particles  in  a  dilute  gas.  This  equation  is  a  powerful  instrument  both  in  the 
analysis  of  thermodynamical  properties  of  fluids  and  in  the  solution  of  fluid  dynamical  problems  in 
rarefied  gas  dynamics,  which  are  often  relevant  in  modern  technology.  The  range  of  problems  of 
rarefied  gas  dynamics  includes,  for  example,  problems  of  flow  past  aircraft  flying  at  high  altitudes, 
motion  of  gases  in  a  vacuum  apparatus,  ultrasonic  vibrations  in  gases,  structure  of  shock  waves, 
nonequilibrium  flows,  etc...  A  survey  of  several  applications  to  typical  fluid  dynamics  flows,  such 
as  shock  waves,  Couette  and  Poiseuille  flows,  etc...  can  be  found  in  the  book  by  Kogan  [36].  The 
ideas  of  kinetic  theory  also  find  applications  in  other  fields,  such  as  radiative  transfer,  the  theory 
of  ionized  gases,  the  theory  of  neutron  transport  and  the  study  of  quantum  effects  in  gases  [10,  20] 

The  Boltzmann  equation  is  a  mathematical  model  of  the  phenomenological  kinetic  theory  of 
gases  which  defines  the  evolution,  in  time  and  space,  of  the  one  particle  distribution  function  for  a 
simple  monoatomic  gas.  The  evolution  equation  for  the  distribution  function 

/  =  /(#,  u,  t)  :  JR3  x  IR3  x  IR+  IR+, 

is  obtained  equating  the  total  derivate  of  /  to  the  gain  and  loss  terms  due  to  the  collisions  of 
the  particles  which  at  the  time  t  €  [0,T]  are  in  a  neighbourhood  of  the  phase  point  x,v,  where 
x  €  D  C  1R3  is  the  position  and  v  €  1R3  is  the  velocity.  The  Boltzmann  equation,  for  a  large  class 


Q(f,f)  (!) 

where  Q+  and  Q-  denote,  respectively,  the  gain  and  loss  terms,  respectively  are  given  by 

Q+(/,  /)  =  /R»xS«  B(z,  ?)/(*,  v,  t)f{x ,  t£,  t)dzdvl 
Q-(f>  f )  =  /(*>  M)  J]R3X52  B(Z,  q)f(z,  v *,  t)dzdv* 

Above 

•  F  is  the  external  force  field  acting  on  the  particles; 

•  z  is  the  unit  vector  in  the  direction  of  the  apse-line  bisecting  velocities  q  =  v*  -  v  and 
/  /  / 

q  =  v*  -  v  ;  ^  b  /  , 

•  v  and  u*  are  the  pre-collision  velocities  of  the  test  and  field  particles,  respectively,  and  v  ,  v-k 
are  the  post-collision  velocities.  Since  for  binary  interactions,  collisions  conserve  ^momentum  and 
energy,  these  velocities  are  related  to  v  and  u*  by  the  relations  v  =  v  +  z(z.q),  «*  =  «*-  2(2.9); 

•  S4.  is  the  integration  domain  of  2,  S+  =  {2  €  1R3  :  \z\  =  1,  <  2, 9  >>  0 1 

9  B  is  a  collision  kernel  which  depends  upon  the  interaction  potential  [20,  16]. 


of  pair-particles  interaction  potentials,  is  [20,  16] 


df 


1  ^  +  vVJ  +  F-V,/ 

at 
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The  Boltzmann  equation  is  a  difficult  model  to  deal  with  both  at  an  analytical  and  computa¬ 
tional  level.  This  difficulty  has  been  an  important  challenge  for  numerous  applied  mathematicians 
also  engaged  by  the  demand  of  aerodynamicists  and  engineers  to  solve  fluid  dynamics  problems 
in  molecular  flow  conditions;  namely  in  physical  conditions  such  that  the  models  of  continuum 
hydrodynamics  cannot  be  any  longer  considered  valid  for  the  applications.  A  great  effort  has  been 
devoted  to  solve  mathematical  and  computational  problems  related  to  the  Boltzmann  equation. 
Some  of  them,  as  we  shall  see,  have  been  (or  partially)  solved  by  several  authors  among  them  the 
author.  These  problems  deal  with  the  following  fundamental  topics 

(i)  the  existence  and  uniqueness  theory, 

(ii)  the  hydrodynamic  limit, 

(iii)  the  boundary  conditions, 

(iv)  the  intermediate  regimes, 

(v)  finding  models  alternative  to  the  one  of  the  classic  Boltzmann  equation. 

Here,  the  intermediate  regimes  refer  to  the  cases  where  the  Knudsen  (Rn)  number  which  mea¬ 
sures  the  ratio  between  the  average  time  separating  two  successive  collisions  of  a  given  particle  and 
a  characteristic  time  of  the  flow  satisfies  either: 

(a)  I<n  <  0.5HT1  or 

(b)  0.510-1  <  Kn  <  10. 

In  the  first  case  it  is  assumed  that  the  Knudsen  number  is  not  extremely  small  ( Kn  «  0.510-1), 
because  otherwise  the  standard  continuous  model  is  valid.  In  both  cases  the  direct  Boltzmann 
simulation  is  very  expensive  and  may  be  impossible  in  case  (a)  because  the  computational  method 
requires  one  cell  per  mean  free  path.  The  hydrodynamic  limit  refers  to  the  study  of  the  asymptotic 
trends  of  the  kinetic  equations  as  the  Knudsen  number  tends  to  zero. 

These  problems  have  been  and  still  (some  of  them)  are  today  a  challenging  objective  and  a  hard 
obstacle  for  applied  and  computational  mathematicians  with  interest  in  this  field  of  mathematical 
physics.  Some  of  them  are  solved  recently  by  the  author,  as  we  shall  see.  A  further  development 
of  his  solution  will  be  considered  in  the  proposed  research. 

Consider  for  instance  the  physically  interesting  problems:  the  interior  domain  problem  which 
is  such  that  the  gas  is  contained  in  a  volume  with  a  solid  surface  on  the  boundary  and  the  exterior 
domain  problem,  which  is  such  that  the  gas  occupies  the  whole  space  1R^  and  contains  an  obstacle. 
The  surface  of  the  solid  wall  is  defined  in  both  cases  by  dTW)  the  normal  to  the  surface  directed 
towards  the  gas  is  v. 

Let  /+  and  f~  denote  the  partial  incoming  and  outgoing  traces  on  the  boundary  dlV  For  / 
continuous,  they  can  be  defined  as  follows 

f  f+{x,v)  =  f{x,v),  xe'drw,  v  ■  n(x)  >  0 
f+(x,v )  =0,  x  €  dTw,  v-n(x)<  0 

and 

f  f~(x,v)  =f(x,v ),  x  <E  dTw,  v  •  n(x)  <  0 
(  f~{x,v)  =0,  x£dVw,  v  ■  n(x)  >  0 

Then  the  boundary  conditions  are  formally  defined  by 

f+(x,v,t)  =  Rf-(x,v,t)  (2) 

where  the  operator  R  is  characterized,  for  a  broad  range  of  physical  problems,  by  the  following 
properties: 
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1)  R  is  linear,  of  local  type  with  respect  to  x  and  is  positive,  i.e.  /  >  0  =>  Rf  (x,  v,  t). 

2)  R  preserves  mass,  i.e.  the  flux  of  the  incoming  particles  equals  the  one  of  the  particles  which 

leave  the  surface. 

3)  R  preserves  local  equilibrium  at  the  boundary,  i.e.  M+  =  RM~ ,  where  Mw  is  the  Maxwellian 
distribution  at  the  wall  temperature. 

4)  R  is  dissipative,  i.e.  satisfies  the  Darrozes  and  Guiraud  [21]  inequality  at  the  wall 

f (vz)(f-  +  Rr)(logf+  ^j-)dv  <  0. 

J  lit 

More  details  about  the  operator  R  and  the  mathematical  aspects  on  the  formulation  of  the 
boundary  conditions  can  be  found  in  [14,  21,  33,  16].  The  experimental  activity  on  the  gas-surface 
interaction  addressed  to  validate  the  mathematical  models  for  the  applications  is  documented  in 
[49,  50].  The  formulation  of  the  initial-boundary  value  problem,  in  the  case  of  the  interior  domain 
problem,  consists  in  linking  the  evolution  equation  (1)  to  the  boundary  conditions  (2)  on  the  wall 
dTw.  In  the  case  of  the  exterior  domain  problems,  in  addition  to  the  boundary  conditions  on  the 
wall,  generally  suitable  Maxwellian  equilibrium  conditions  at  infinity  are  assumed. 

Mathematical  problems  related  to  the  applications  may  involve  additional  difficulties.  For 
instance  it  may  be  necessary  to  deal  with  non  regular  or  concave  geometries.  In  addition,  it  may  be 
useful  or  even  necessary,  to  decompose  the  domain  occupied  by  the  gas  into  domains  where  one  can 
apply  the  equations  of  the  continuum  fluids  and  into  domains  where  one  has  to  apply  the  Boltzmann 
equation.  Such  decomposition  has  been  introduced  by  the  author  [73].  It  has  several  computational 
advantages,  as  we  shall  see,  but  involves  additional  mathematical  and  physical  difficulties  related 
to  the  matching  of  the  two  models. 

Since  the  existence  theory  involves  rather  restrictive  assumptions,  it  is  only  partially  useful  for 
the  applications.  The  most  important  existence  results  are  obtained  by  DiPerna  and  Lions  [23,  24] 
who  considered  the  solutions  for  large  L1  initial  data.  This  type  of  results  exploit  the  averaging 
Lemma  by  Golse,  Lions,  Perthame  and  Sentis  [28]  to  obtain  existence  of  weak  renormalized  solu¬ 
tions.  A  further  important  developments  are  obtained  in  [34,  2,  37,  40].  For  a  complete  survey  on 
the  recent  development  on  the  existence  theory  the  reader  is  refered  to  [3]. 

The  other  topics  (ii)-(v)  were  (and  are  still)  covered  to  some  extent  by  the  author  [73]-[88],[llJ- 
[12]  ,[51]-[56].  The  program  proposed  here  is  a  continuation  of  the  prior  and  current  research  of 
the  author.  Since  it' will  use  mainly  the  methods  developed  by  the  author,  his  prior  and  current 
research  are  included  with  a  little  more  detail  and  the  rest  of  this  proposal  will  refer  mainly  to  the 
work  by  the  author  [73]-[88],  [ll]-[12],[51]-[56]. 

2  Current  and  prior  research  of  the  author  related  to  (ii)-(v) 

The  Boltzmann  equation  is  one  of  the  most  important  tools  in  gas  dynamics  calculations,  when 
physical  phenomena  of  a  molecular  scale  cannot  be  neglected.  In  this  case  the  model  of  continuum 
hydrodynamics  cannot  be  any-longer  considered  valid  for  the  applications.  But,  when  the  mean  free 
path  gets  too  small,  the  numerical  solution  of  Boltzmann  equations  becomes  impossible  because  the 
discretization  step  of  the  associated  grids  must  be  smaller  than  the  mean  free  path.  The  classical 
solution  consists  then  in  replacing  the  Boltzmann  equation  by  its  fluid  limit  obtained  when  the 

mean  free  path  goes  to  zero.  . 

The  hydrodynamical  limit  theory  aims  to  find  the  connection  between  the  Boltzmann  equation 

and  the  hydrodynamical  equations  (Euler  and  Navier-Stokes  equations).  This  can  be  formulate 
as  the  search  of  asymptotic  relationships  between  solutions  of  these  equations.  This  connection 
results  from  two  types  of  properties  of  the  collision  operator: 
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1)  Conservation  properties  and  an  entropy  relation  that  implies  that  the  equilibrium  distribution 
corresponds  to  a  Maxwellian  for  the  limit  at  the  order  zero. 

2)  The  derivative  of  the  collision  operator  satisfies  a  formal  Fredholm  alternative  with  a  kernel 
related  to  the  conservation  properties  of  1). 

The  compressible  Euler  equations  are  obtained  formally  using  the  conservation  properties  and 
the  entropy  dissipation  which  are  consequences  of  the  properties  of  the  collision  operator.  In  the 
Chapman-Enskog  or  Hilbert  expansion  [20]  of  fin  €  =  I<n  the  compressible  Euler  equations  are  the 
leading-order  dynamics. 

Since  the  compressible  Euler  equations  generally  become  singular  after  a  finite  time  [71],  any 
global  in  time  convergence  proof  cannot  rely  on  uniform  regularity  estimates.  In  [7]  assumptions 
are  made  on  the  kinetic  level.  The  authors  assume  a  formally  consistent  convergence  for  the  fluid 
dynamical  moments  and  entropy  of  the  solution  of  the  kinetic  equations  (1).  A  more  detailed 
knowledge  of  the  collision  operator  is  needed  in  order  to  obtain  the  compressible  Navier-Stokes 
equations.  These  equations  arise  as  corrections  to  those  of  Euler  at  the  next  order  in  the  Chapman- 
Enskog  expansion.  Strong  hypothesis  are  needed  on  the  regularity  of  solutions  of  the  compressible 
Navier-Stokes  equations  in  order  to  make  sens  of  these  expressions.  The  results  available  up  to  now 
consider  only  the  case  of  the  full  space  or  periodic  domain  [7,  8,  15,  60]. 

In  presence  of  obstacles  boundary  conditions  must  be  specified.  When  the  Knudsen  number 
is  extremely  small  (very  dense  gas),  the  classical  hypothesis  of  no-slip  boundary  conditions  give 
accurate  boundary  conditions.  However,  when  the  Knudsen  number  lies  in  the  range  of  values 
of  the  case  (a)  and  (b)  (iv),  a  breakdown  of  the  aerodynamic  theory  in  the  region  neighboring 
the  obstacle  (the  so-called  Knudsen  layer)  is  observed  and  boundary  conditions  of  slip  type  must 
be  specified.  The  Chapman-Enskog  expansion  (or  any  other  method  based  on  the  use  of  a  finite 
number  of  moments)  does  not  in  general  satisfy  the  kinetic  boundary  conditions.  Therefore,  this 
expansion  is  not  valid  in  the  Knudsen  layer  which  is  of  width  of  the  order  oi  the  mean  free  path. 

The  standard  solution  is  to  use  analytical  slip  boundary  conditions  as  described  in  [36,  13,  16]. 
But  the  constants  which  are  involved  are  hard  to  identify  and  their  validity  is  questionable.  On 
the  other  hand,  the  direct  simulation  of  the  kinetic  problem  is  rapidly  too  expensive,  because  it 
requires  one  computational  cell  per  mean  free  path.  To  overcome  such  difficulties  many  authors 
have  recently  tried  to  use  intermediate  asymptotic  models  such  as  Burnett  equations  [89].  However, 
in  [18]  the  authors  have  shown  that  the  Burnett  equations  are  in  violation  of  the  second  law  of 
thermodynamics  thus  explaining  their  long  history  of  numerical  difficulties. 

On  the  other  hand  when  the  mean  free  path  is  roughly  one  thousand  times  smaller  than  the 
length  of  the  obstacle  these  asymptotic  models  are  no-longer  valid. 

In  [73],  the  author  proposed  a  fundamental  strategy  that  permits  the  coupling  of  differents 
models  and/or  different  approximations  to  compute  the  solution  of  the  exterior  domain  problem. 
This  strategy  is  very  general  and  can  be  applied  to  many  different  physical  systems  [73]-[88],  [51]- 
[56],[ll]-[12].  The  application  of  this  strategy  to  the  solution  of  intermediate  regimes  (as  defined 
in  (iv))  consists  of  coupling  the  hydrodynamics  equations  (Euler  or  Navier-Stokes  equations)  with 
the  Boltzmann  equation.  The  resulting  method  involves  additional  mathematical  and  physical 
difficulties  related  to  the  matching  of  equations  of  the  two  models.  However,  this  approach  has 
several  computational  advantages.  One  of  the  great  advantages  is  the  use  of  the  correct  model 
related  to  the  physical  features  of  the  flow. 

The  application  of  this  approach  to  the  solution  of  Boltzmann  equation  for  external  domain 
problem  consists  of  the  following  steps 

(1)  Domain  decomposition  of  the  external  field  into  the  domains,  possibly  overlapped,  of  validity 
of  the  Boltzmann  Equation  and  of  the  equations  of  hydrodynamics. 

(2)  Solution  of  the  kinetic  equations  in  the  domain  of  validity  of  the  Boltzmann  equation. 
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(3)  Solution  of  the  hydrodynamics  equations  in  their  domain  of  validity. 

(4)  Coupling  of  the  solution  to  the  two  models,  the  continuous  and  kinetic  ones. 

These  steps  are  documented  in  the  next  paragraphs.  The  first  paragraph  deals  with  the  asymp¬ 
totic  models  at  the  fluid  limit.  The  description  of  the  geometry  and  the  coupling  of  these  asymptotic 
models  with  the  original  Boltzmann  equation  is  then  discussed  in  the  second  paragraph.  The  third 
paragraph  presents  a  brief  description  of  the  different  discretization  methods  employed  for  the 
kinetic  model  and  the  hydrodynamic  one.  In  the  fast  paragraph  some  comments  are  presented. 

2.1  Asymptotic  models 

The  Boltzmann  equation,  in  dimensionless  form,  can  be  written  as  follows  [16] 

%  +  v-Vxf=-Q(f,f).  (3) 

at  e 

The  number  e  =  Kn  is  the  Knudsen  number  defined  in  the  previous  section.  The  fluid  regime 
corresponds  to  very  small  values  of  e  (e  <  0.01).  In  such  regimes,  the  governing  equations  can  be 
obtained  by  a  formal  expansion  of  Boltzmann  equation  in  power  of  e.  At  first  order,  the  Chapman- 
Enskog  expansion  assumes  that  /  is  as  follows 


/  =  /o(l  +  £<t>)- 

(4) 

By  writing  the  Boltzmann  equation  at  first-order  in  e  yields 

-eQ(foJo)  =  0. 

(5) 

The  distribution  /0  is  therefore  a  collision  invariant,  that  is  a  Maxwellian 

fo(x,v,t)  =  nu€(u{x,t)). 

(6) 

Plugging  this  value  of  fo  into  the  Boltzmann  equation  yields 

-Jr  +  v-  Vj;/  =  Q(fo ,  /o<£)  +  Q(fo<f>,  fo)  +  0{6)’ 

^  ut 

and  introducing  the  linear  operator  L<f>  =  Q(/o>  /o^)  d"  /o)?  yields 

(7) 

II 

*|§ 

+ 

<1 

H 

0* 

II 

+ 

Q 

n\ 

(8) 

that  is  (assuming  L  to  be  invertible) 

(9) 

The  dependence  of  /0  in  x  is  obtained  by  solving  the  conservation  laws,  that  follows  from  the 
Boltzmann  equation  multiplied  by  U  =  (M,  M2)  and  integrated  in  the  velocity  domain.  This 
integration  cancels  the  collision  terms  (since  mass,  momentum  and  energy  are  conserved  during 
collision)  and  yields 

Ft  Ju3  /o(1  +  €®udv  +  VxIjr3  v/o(1  +  €<t>)Udv  =  °'  (10) 

Up  to  now,  all  equations  are  exact.  Three  simplifications  are  then  made  in  the  explicit  calcula¬ 
tion  of  4>.  In  particular 
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i)  The  left  hand  side  term  in  Eq.  (9)  is  computed  at  first  order  by  writing  the  conservation  laws 
(10)  at  order  0  in  6. 

ii)  The  correction  term  <f>  is  normalized  by  means  of  the  condition 


In  other  words,  the  parameters  p)  u  and  T  of  the  Maxwellian  /o  correspond  to  the  fluid  macro¬ 
scopic  density,  velocity  and  temperature, 

iii)  The  perturbation  equation 


Ld>=^  +  0(e)  (12) 

at 

is  only  solved  within  the  subspace  generated  by  the  first  Laguerre-Sonine  polynomials. 

With  these  assumptions,  the  approximate  solution  of  the  Boltzmann  equation  (written  at  first 
order  level)  becomes 


/  —  fcE 


where 


n  ,  crNM  *  *  /lcl  V7/3 

(27rr^)feXP  "  2  )[1  5p{r6)^2re  2)C'  p{r6)^ 

t  |c|2  2T-c  c2  5  1A  ;c®c 

(27rr0)ieXP(  2)[1+5p(r9)*{2r6  2>  2  p(r0?  J’ 


-c2Id)  :  Vuj 


A  =  —  j  mfoccfrc  ®  cdc  — 

p[V«  +  V0  —  Wuld], 

(14) 

T  =  J  foe4>mc2cdc  = 

-xve, 

(15) 

c  = 

V  —  u 

(16) 

Vrd  ’ 

A  and  p.  are  viscosity  coefficients  depending  of  the  form  on  the  collision  kernel  and  on  the  degree 
of  the  Sonine  plynomials  used  in  the  calculation  of<p,  and  Id  denotes  the  indentity  tensor  and  A  :  H 
the  quantity  Tr(AH). 

Plugging  this  expression  of  f  back  in  the  conservation  laws,  finally  gives  the  Navier-Stokes 
equations 

p  1  T  pu  1 

pu  +V*  pu®u  +  pld  =0  (17) 

pE  ((pE  +  p)Id-A).u  +  q 

The  conservation  laws  Eq.  (10)  are  quite  natural  from  the  physical  point  of  view  and  are 
justified  on  the  mathematical  level.  Indeed,  by  writing  the  Boltzmann  equation  for 

/  =  fcE  +  **  h  —  /o(l  +  €0)  +  f2/2!  (1®) 

one  gets,  at  order  2  in  e 
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dfcE 

dt 


=  Q(/o>  M  +  e/2)  +  <3(/o</>  +  e/2,/o)  +  eQ(M,  /o</>) 


=  L((j>  +  +  tQifofi)  fofy  * 

JO 

(19) 

This  equation  has  a  solution  of  the  type  <^+  if  and  only  if  the  term 

JO 

dfdt  eQ^M) 

(20) 

belongs  to  the  image  of  L ,  i.e. 

(df°tE  eQtfohfoM)  -  0,  (ImZ/)T, 

(21) 

with  (Im£)T  of  the  form 

(ImI)T  =  {v>,  V>  =  a  +  /3u  +  7^M2} 

(22) 

Because  of  the  invariance  properties  of  the  collision  operator  Q,  this  solvability  condition  reduces 
exactly  to  the  conservation  laws  (10). 


2.2  Coupling 

Consider  the  flow  of  rarefied  gas  in  a  global  domain  around  a  solid  body  of  boundary  F w ■  This 
domain  is  splitt  into  two  (possibly  overlapped)  subdomains  and  The  domain  Qb  will  be  used 
for  the  Boltzmann  equations.  This  is  a  local  domain  which  contains  the  obstacle,  with  internal 
boundary  and  external  boundary  T t,  C  &e-  The  domain  Qe  will  t*e  use(f  f°r  tl*e  Euler  or  Navier- 
Stokes  equations.  This  is  a  large  domain,  with  an  internal  boundary  r,n<  C  which  surrounds  the 
body  and  an  external  boundary  reit  which  is  the  external  boundary  of  the  computational  domain 
Q.  Depending  of  the  gas  density,  either  the  Boltzmann  domain  Qb  or  the  fluid  domain  Qe  can  be 

suppressed. 

We  shall  now  describe  a  numerical  method  for  the  simultaneous  solution  of  the  Navier-Stokes 
and  Boltzmann  equations  when  these  equations  are  used  together  but  on  two  different  domains. 
For  this  purpose,  we  first  introduce  compatibility  conditions  to  be  imposed  at  the  interfaces,  and 
we  then  propose  an  algorithm  for  the  numerical  solution  of  the  resulting  coupled  problem. 

In  a  first  case  where  F,nt  =  F„  (Qb  C  fte:  full  overlapping),  we  must  specify  boundary  conditions 
for  the  Navier-Stokes  equations  on  the  wall  Tu,.  These  boundary  conditions  cannot  be  obtained 
from  the  asymptotic  model  since  the  latter  is  based  on  the  Chapman-Enskog  model  which  is  not 
valid  in  the  kinetic  layer.  These  are  obtained  by  considering  boundary  fluxes  calculated  using  their 
discrete  kinetic  definition 


"  —  dVKAt 

■  v-  ~ 


ZieM  ~  »f ) 


(23) 


which  correspond,  respectively  to  momentum  and  energy  fluxes.  Here,  /  denotes  the  set  of  particles 
which  collide  with  the  subset  dr  of  the  wall  during  the  time  KAt,  and  v{  (respectively  vf)  denotes 
the  pre- (respectively  post-)  collision  velocity  of  the  individual  particle  i  of  mass  m. 

In  the  second  case  where  rint  #  Tw  the  compatibility  conditions  are  simpler  to  obtain  and  are 
based  on  the  above  kinetic  derivation  of  the  Navier-Stokes  equations.  Therefore,  on  each  domain, 


8 


we  only  need  to  impose  as  boundary  conditions  the  kinetic  distribution  of  the  incoming  particles. 
For  the  fluid  domain,  this  ditribution  is  taken  to  be  equal  to  the  distribution  fB  of  particles  leaving 
the  Boltzmann  domain  at  the  same  point.  Similarly,  for  the  Boltzmann  domain,  the  distribution 
of  incoming  particles  is  taken  to  be  equal  to  the  distribution  of  particles  fcE  leaving  the  fluid 
domain  at  this  point,  where  fcE  Is  computed  from  the  Navier-Stokes  macroscopic  quantities,  by 
the  relation  (13).  Thus,  in  the  Boltzmann  domain  tlb,  the  Boltzmann  equation  (1)  is  solved  with 
the  boundary  conditions  (2)  at  the  wall  and 

f{x,v,t)  =  fcE{x,v,t) ,  \/vuB<  0,  (24) 

on  the  Navier-Stokes  interface  Tb.  In  the  fluid  domain  £2ei  the  Navier-Stokes  equations  (17)  are 
solved.  In  these  equations,  the  kinetic  boundary  condition  fcE{x,v,t)  =  fB(x,v,t),  Vu  •  ue  <  0 
cannot  be  imposed  pointwise  and  must  therefore  be  integrated  in  v,  yielding  the  Neumann  type  of 
boundary  condition 


pu 

pu®u-\-  pld  —  A 
{{pE  +  p)u  +  q  -  A  :  Vu 


v 


F+  +  F~ 


(25) 


1 

r 

1 

V 

fcE  + 

/  ( V-v ) 

V 

.  IM2 . 

J  V'i/<  0 

.  §M2 . 

The  coupling  algorithm  is  the  transmission  time  marching  algorithm  introduced  and  studied  by 
the  author  in  [73].  It  successively  solves  the  Boltzmann  equations  on  the  Boltzmann  domain  with 
imposed  boundary  conditions  deduced  from  the  previous  solution,  and  the  Navier-Stokes  equations 
on  the  fluid  domain  with  imposed  boundary  conditions  deduced  from  the  Boltzmann  solution. 

For  n  >  0,  and  knowing  estimates  (pn,pun,En)  and  /£  of  the  Navier-Stokes  and  Boltzmann 
solutions  on  Qe  and  new  estimates  /£+1  and  {pn+1 ,  pun+1 ,  En+l)  are  therefore  computed  by: 

i)  solving  the  Boltzmann  equation  (1)  on  Q,b  with  initial  condition  /(.,0)  =  fB  and  boundary 
conditions  (2).  At  the  end  of  this  step,  we  obtain  an  improved  estimate  fg+1  of  the  Boltzmann 
solution  on  £2e  and  an  updated  value 


fg+1(x,v,e)dv  (26) 

of  the  half-flux  entering  through  r,n<; 

ii)  solving  The  Navier-Stokes  equations  (17)  on  Qe  with  boundary  conditions  (25),  where  the 
value  of  the  half-flux  entering  Qe  through  rinf  is  replaced  by  the  value  (Fn+1)“  computed  in  Eq. 
(26)  if  Tint  7*-  r^,,  and  boundary  conditions  (23)  if  Tint  = 

In  the  sequel  the  coupling  with  Tint  =  Tw  will  be  refereed  to  as  the  coupling  (a)  and  the 
coupling  with  T^t  7^  ru,  will  be  refereed  to  as  the  coupling  (/?)• 


(Fn+ir  =  f  ( 

Jvv<0 


{vu) 


L IM2 


2.3  Discretization 

For  the  Boltzmann  domain  the  discretization  method  is  based  on  splitting  the  physical  processes 
into  a  free  transport  phase  and  a  collision  phase.  This  splitting  converges,  as  the  time  step  tends 
to  zero  [22]  to  the  DiPerna-Lions  [23,  24]  renormalized  solution  of  the  Boltzmann  equation.  The 
collision  phase  is  solved  by  so-called  random  particle  methods  which  are  more  closely  related  to 
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the  Boltzmann  equation  [4]-[6],[46]-[48].  The  analysis  of  convergence  of  the  particle  method  can 
be  found  in  [4]-[6],[66].  For  the  fluid  domain  three  methods  have  been  used:  the  hybrid  finite 
element/finite  volume  methods  [68,  79],  the  streamline  upwind  Petrov-Galerkin  methods  [44,  69] 
and  the  hybrid  upwind  splitting  [19].  The  coupled  approach  introduced  new  boundary  conditions 
and  geometry  and  their  treatment  constitutes  an  important  step  in  the  discretization  process  [73, 
12].  -  - 

2.4  Comments 

Compared  with  the  coupling  strategy  (/3)  which  uses  very  little  overlapping  between  the  kinetic 

and  hydrodynamic  region,  the  coupling  (a)  is: 

-simpler  because  the  Navier-Stokes  equations  are  used  in  the  whole  physical  domain, 

-less  sensitive  to  the  choice  of  the  interface  rtn<,  which  makes  it  more  robust  at  dense  regimes 
especially  in  the  case  where  the  direct  Boltzmann  simulation  is  not  possible  because  it  requires  one 
computational  cell  per  mean  free  path. 

-but  restricted  to  rather  dense  situations  (case  (a)  of  Section  1)  because  it  integrates  the  hy¬ 
drodynamic  model  up  to  the  wall. 

-gives,  without  changing  the  global  Navier-Stokes  solver,  an  easy  way  of  supplementing  and 
testing  a  large  variety  of  kinetic  boundary  conditions.  These  kinetic  boundary  conditions  are  first 
imposed  on  the  Boltzmann  model,  and  the  resulting  fluxes  crB  and  qB  are  then  plugged  in  the 
Navier-Stokes  equations.  They  correspond  to  the  losses  in  tangential  momentum  and  energy  of 

particles  colliding  into  the  wall  [76].  _  ... 

By  contrast,  the  coupled  approach  (ft)  is  used  mainly  for  rarefied  regimes,  and  faces  difficulties 

for  denser  cases. 

Therefore,  the  solution  procedure  described  herein  provides  a  fundamental  contribution  to  the 
solution  of  the  problem  of  intermediate  regimes  (iv),  the  problem  of  boundary  conditions  for 
hydrodynamic  limits  (iii)  and  to  some  extent  the  problem  (v).  Moreover,  this  solution  procedure 
has  been  successfully  developed  not  only  for  a  simple  gas,  but  also  for  particles  with  internal  degrees 
of  freedom.  In  particular,  the  coupling  strategy  allowed  the  treatment  of  experimental  situations  in 
the  intermediate  regimes.  The  results  obtained  are  compatibles  with  available  experimental  results 
[51]-[52],[12].  Hence,  they  provide  a  major  contribution  to  the  solution  of  transitional  regimes. 
More  details  about  t\ie  different  aspects  of  the  solution  procedure  presented  here,  including  its 
mathematical,  numerical  and  physical  validity,  can  be  found  in  the  series  of  papers  by  the  author 

[73]-[88],[51]-[56],[ll]-[12].  .  ,  ,  .  4l  . 

Although  the  solution  procedure  described  above  is  very  flexible,  the  author  is  currently  inves 
tigating  the  possibility  of  making  it  even  more  flexible.  For  instance  the  determination  of  a  suitable 
simple  criterion  in  order  to  choose  an  optimal  Boltzmann  domain.  Also,  certain  difficulties  met  in 
the  study  of  these  methods  need  to  be  removed.  For  instance,  if  we  compare  the  CPU  time  of  the 
coupling  approach  and  the  asymptotic  model  (in  the  range  of  its  validity:  Kn  <10  4),  the  coupling 
strategy  requires  more  than  twice  the  CPU  time  of  a  direct  asymptotic  model  simulation.  However, 
the  latter  (within  its  range  of  validity: I<n  <  10“4)  fails  for  more  complex  physical  situations,  while 
the  coupled  strategy  works  well.  Also  the  Navier-Stokes  calculations  performed  within  the  cou¬ 
pling  approach  required  robust  solvers  and  adequate  meshes.  This  led  the  author  to  pursue  further 
investigation  of  the  solution  procedure.  For  instance,  the  developement  of  new  coupling  models  is 
considered.  This  has  to  be  done  in  connection  with  the  problem  (v).  In  the  solution  procedures 
developed  above  one  has  to  deal  with  a  singular  perturbed  problem,  where  the  small  parameter 
is  the  mean  free  path.  In  [7,  8],  the  different  incompressible  limits  of  the  Boltzmann  equations 
are  obtained  using  a  very  interesting  scaling  of  both  the  initial  data  and  the  equation.  The  small 
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parameter  they  use  is  the  ratio  between  the  Mach  number  and  the  Reynolds  number.  It  would  be 
interesting  to  develop  the  coupling  techniques  according  to  different  scaling.  Finally,  there  is  no 
theoretical  estimate  of  the  number  of  time  steps  within  each  model  calculation  for  each  coupling 
iteration.  This  motivates  further  mathematical  development  related  to  these  coupling  methods. 
Some  of  these  constitute  part  of  the  proposed  research  as  documented  in  the  next  sections. 

3  Development  of  new  coupling  models 

The  objective  here  is  to  develop  new  coupling  models.  From  the  previous  section  it  is  clear  that 
the  most  difficult  part  to  solve  in  the  intermediate  regime  corresponds  to  case  (a)  of  Section  1. 
The  direct  Boltzmann  simulation  becomes  impossible  near  the  limit  range  of  validity  of  the  direct 
asymptotic  approach.  The  solution  procedure  described  above  provides  an  efficient  solution  to  this 
problem,  but,  becomes  expensive  in  the  limit  range  of  validity  of  the  direct  asymptotic  approach. 
The  objective  is  then  to  find  a  model  that  would  replace  the  Boltzmann  equation.  This  model  should 
retain  the  pertinent  physical  properties  of  the  gas  while  having  simpler  mathematical  form.  The 
BGK  model  introduced  by  Bhatnager,  Gross  and  Krook  [10]  satisfies  these  requirements.  Hence, 
this  is  the  model  that  will  be  considered  here.  The  main  difficulty  in  handling  the  Boltzmann 
equations  arises  from  the  complicated  nature  of  the  collision  terms.  The  BGK  model  consists 
precisely  in  replacing  the  collision  integrals  by  mathematically  simpler  terms;  these  however,  are 
chosen  so  as  to  conform  to  the  conservation  laws  for  mass,  momentum,  and  energy  and  at  the 
same  time  to  represent  certain  essential  features  of  collisions.  It  corresponds  to  the  following  set  of 
equations 

%  +  „.?,/  =£(/„-/)•  (27) 

at  \x 

where  P  is  the  pressure,  p  is  the  Navier-Stokes  viscosity  and 

p  |u  -  ft]2 

f°  =  (2tt i?T)3/2eXP(  2RT  ’ 

is  the  local  Maxwellian.  The  macroscopic  quantities  are  related  to  /  through  the  moments: 

p{x,  t )  =<  /  >,  u(z,  t)  =  i  <  vf  >,  RT=  <  lcl 2f  >> 

c  —  v  u,  <  /  >=  J^3f{x,t,v)dv. 

The  appearance  of  the  local  Maxwellian  /0  means  that  the  equation  remains  a  nonlinear  integro- 
differential  equation.  The  existence  of  a  solution  to  this  nonlinear  equation  has  been  established  by 
Perthame[62].  This  approach  makes  possible  a  survey  of  the  whole  range  from  low  density  to  high 
density,  including  the  intermediate  regime.  Moreover,  it  satisfies  the  correct  microscopic  boundary 
conditions.  This  is  fundamental  for  the  objective  of  the  actual  proposed  research.  The  model  leads 
also  to  the  correct  asymptotic  behavior  in  the  two  limiting  cases:  Kn  — t  0  and  Kn  — f  oo.  Moreover, 
this  method  is  capable  of  considerable  extension  and  generalization  to  include  physically  complex 
problems. 

This  model  has  the  advantage  of  describing  the  right  fluid  limit.  However,  the  transport  coeffi¬ 
cients  obtained  at  the  Navier-Stokes  level  through  Chapman-Enskog  expansion  are  not  satisfactory. 
For  instance  the  Prandtl  number  Pr  =  relating  the  viscosity  v  to  the  heat  conduction  k 
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is  equal  to  1.  For  most  gases  Pr  <  1.  Using  Boltzmann  equation  the  hard-sphere  model  for  a 
monoatomic  gas  yields  a  Prandtl  number  =  2/3.  Several  modification  to  the  BGK  model  have 
been  proposed  [13,  16].  These  however  do  not  satisfy  simultaneously  the  requirements  of  yielding 
a  nonegative  distributions,  predicting  a  Prandtl  number  less  than  1,  and  satisfying  the  celebrated 
Boltzmann  theorem:  f  Q(f) In fdv  <  0.  A  model  satisfying  the  first  two  requirements  was  proposed 
by  Holway  [31].  The  third  requirement  was  proved  to  be  true  for  the  Holway  model,  also  known 
as  the  Gaussian  or  ellipsoidal  (ES-BGK)  model,  only  very  recently  [1].  The  Gaussian  BGK  model 
corresponds  to 


where 


dj_ 

dt 


+  v-Vxf 


P 

H(1  -  u) 


(. fo~f )• 


(28) 


fo  = 


\/det(27rA) 


exp(- 


--(u-  u)  •  A 


-l 


(u-  u)), 


A  =  (1  -  v)RTM  +  uQ, 
T=4;tr0>  0  =  -  <  c<g>c/>, 


and  the  parameter  <  v  <  1  is  introduced  to  give  the  right  Prandtl  number  through  the 
formula:  |  <  Pr  =  —  +°°- 

The  new  proposed  coupling  model  is  therefore  to  use  the  Gaussian  BGK  (ES-BGK)  .model  in 
place  of  the  Boltzmann  model  in  the  coupling  strategy  introduced  in  the  previous  sections.  Since 
the  Chapman-Enskog  method  may  be  applied  to  the  ES-BGK  model,  the  compatibility  conditions 
and  the  coupling  (a)  and  (/?)  for  the  Gaussian  BGK  are  derived  in  a  similar  fashion  as  in  Section 
2.  The  study  of  the  resulting  algorithms  includes  the  analysis  of  the  coupling  algorithms,  and  the 
numerical  study  of  the  solution  procedures. 


3.1  Numerical  studies  of  the  coupling  (a)  and  ( [3 )  for  the  Gaussian  BGK 

A  first  step  in  the  numerical  study  of  our  strategies  for  the  Gaussian  BGK  model  consists  in 
validating  this  model  by  performing  comparison  tests  with  experiments  and  against  the  results 
obtained  using  the  Boltzmann  model.  To  do  this,  we  first  used  a  variant  of  the  Direct  Simulation 
Monte  Carlo  method,  which  we  have  already  developed  and  studied  for  the  full  Boltzmann  model 
[12].  We  have  also  studied  the  Gaussian  BGK  model  for  polyatomic  gases.  The  detail  of  this  study 
is  given  in  [88].  As  an  illustration  of  this  we  present  the  numerical  results  for  the  following  test 
cases.  The  first  case  corresponds  to  a  computation  of  a  flow  around  an  ellipse  for  a  semi-rarefied 
gas  with  the  following  data 


Monoatomic  case  with  R  =  287 
Mean  free  path  Aoo  =  8.310-4 

Characteristic  length  L  =  0.1m,  which  yields  Kn^  =  8.310 
Poo  =  5.1910-5fcp/m3 
^  =  14.5K 

Voo  =  1663 m/s,  which  yields  Machoo  =  20 

At  the  wall  Tw  =  290 K  and  we  use  total  accomodation 
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Viscosity  law:  /j,  =  Moo(^“)a<  with  w  =  0.657  and  =  0.24910-5Pa.s.  At  the  wall  = 
0.17810_4Pa.s.  This  yields  Re oo  =  3469  and  Rew  =  484. 

Computational  domain:  0.12m  x  0.08m. 

Rectangular  mesh:  hx  =  10-3m,  hy  =  2.510-4to 

Time  step:  At  =  0.20410-6s 

Number  of  particles  in  a  cell  at  infinity:  25 


The  numerical  results  for  this  case  are  presented  in  Figures  1-5.  The  second  test  case  corresponds 
to  a  flow  around  an  ellipse  in  the  transitional  regime  case.  This  corresponds  to 

Mean  free  path  Aqq  =  1.010-4 

Characteristic  length  L  —  0.1m,  which  yields  Kn^  =  1.010-3. 

Computational  domain:  0.12m  X  0.08m. 

Rectangular  mesh:  hx  =  8.010-4m,  hy  =  2.010_4m 
Time  step:  At  =  0.1610-6s 


The  remaining  data  are  the  same  as  in  the  first  test  case.  The  numerical  results  for  this  case 
are  presented  in  Figures  6-10.  The  third  test  case  corresponds  to  a  flow  around  an  ellipse  in  the 
transitional  regime  case  for  a  polyatomic  gas.  The  data  are  the  same  as  in  the  second  test  case. 
The  numerical  results  are  shown  in  Figures  11-15. 

These  results  together  with  our  study  [88]  show  that  in  the  transitional  regime  the  Gaussian 
BGK  model  gives  results  that  are  in  good  agreements  with  those  given  by  the  Boltzmann  model. 
This  clearly  justifies  the  study  of  our  coupling  strategies  in  which  the  Boltzmann  model  is  replaced 
by  the  Gaussian  BGK  model.  We  are  currently  performing  a  numerical  study  of  our  coupling  (a) 
and  (ft)  for  the  Gaussian  BGK  model.  In  these  coupling,  we  are  using  DSMC  method  for  the 
discretisation  of  the  Gaussian  BGK,  while  for  the  Navier  Stokes  computation  we  are  using  the 
same  numerical  methods  used  already  for  the  discretisation  of  the  Navier-Stokes  equations  in  the 
coupling  (<*)  and  (/?)  for  Boltzmann  model  described  in  section  2.3.  We  are  considering  several 
test  cases  where  experimental  data  are  available.  We  shall  also  apply  our  methods  to  experimental 
situations  where  the  direct  and  classical  methods  are  not  feasible.  This  will  show  the  fundamental 
importance  and  contribution  of  our  strategies. 

In  order  to  simplify  further  the  numerical  procedure  and  cut  down  the  cost  in  terms  of  com¬ 
puter  time  and  memory,  we  are  investigating  deterministic  approaches  for  the  Gaussian  BGK.  The 
polyatomic  Gaussian  BGK  model  was  obtained  using  an  internal  energy  parameter  /  which  takes 
into  account  the  degrees  of  freedom  in  a  general  way  just  like  the  Boltzmann  case  [12]  (and  the 
references  therein).  Numerically  this  is  expensive.  For  BGK  models,  the  introduction  of  an  internal 
energy  variable  is  not  necessary.  In  fact  all  macroscopic  variables  can  be  equivalently  obtained  by 
a  simpler  procedure.  It  consists  of  introducing  two  density  distributions,  one  for  mass  and  one  for 
internal  energy.  We  obtain  a  BGK  system  governing  the  evolution  of  the  two  distributions  [88]  (and 
the  references  therein).  We  shall  perform  a  thorough  numerical  study  for  this  new  BGK  system. 
We  shall  first  extend  the  Direct  Simulation  Monte  Carlo  (DSMC)  method  mentioned  above  (  for 
the  case  of  polyatomic  gase  using  an  internal  energy  variable).  We  shall  then  compare  the  results 
to  a  DSMC  computation  for  the  (polyatomic)  Gaussian  BGK  using  an  internal  energy  variable. 
We  shall  also  extend  the  deterministic  approach  we  are  intending  to  develop  for  the  monoatomic 
Gaussian  BGK  model  to  the  polyatomic  case.  Another  numerical  investigation  that  we  are  planning 
to  accomplish  is  the  numerical  study  of  the  Gaussian  BGK  for  mixtures. 
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3.2  Theoretical  study  of  the  methodologies  (a)  and  (/?) 

To  establish  the  mathematical  foundations  of  the  proposed  models  (a)  and  (/?),  we  shall  look,  in 
particular,  at  the  following  topics: 

(i)  Hydrodynamic  limits  and  the  rigorous  derivation  of  the  models  (a)  and  (/?). 

(ii)  Existence,  stability,  regularity,  and  positivity  of  the  solutions  of  the  models  (a)  and  (/3). 

(iii)  Asymptotic  analysis  and  trend  to  equilibrium  for  the  models  (a)  and  ((3). 

(iv)  Convergence  analysis  of  the  approximations  of  the  models  (a)  and  (/3)  by  the  transmission 
time  marching  algorithm  and  error  estimates. 

We  shall  present  below  our  recent  accomplishments  together  with  our  future  plans  concerning 
each  of  the  above  points  (i)-(iv). 

(i) :  The  hydrody'namical  limit  is  of  fundamental  importance;  it  is  at  the  heart  of  the  develop¬ 
ment  of  our  models.  The  foundations  of  the  hydrodynamic  limit  has  been  established  by  Hilbert 
for  kinetic  models  of  rarefied  gases  in  his  famous  paper  [35].  This  method  has  been  later  improved 
by  Chapman  [17]  and  Enskog  [26]  to  account  for  dissipative  terms.  A  mathematical  framework  has 
been  provided  by  Grad  [29].  Most  of  the  existing  results  in  this  direction  are  obtained  for  simple 
models  and  periodic  or  full  space  problems.  In  our  case,  we  want  to  solve  problems  in  general 
domains  and  with  general  boundary  conditions.  The  first  simplified  Boltzmann  model  we  studied 
is  a  Carleman  model.  For  this  model  we  proved  the  global  convergence  of  Chapman-Enskog  on  a 
bounded  domain  [80].  This  is  the  first  result  of  hydrodynamic  limit  on  bounded  domain.  Using 
this  result  we  have  been  able  to  derive  rigorously  the  coupling  (ex)  and  (/?)  for  the  Carleman  model. 
We  have  proved  similar  results  for  the  linearized  Boltzmann  model  [81].  In  a  very  recent  work  we 
have  performed  a  rigorous  derivation  of  the  coupling  (cx)  and  (/?)  for  a  (nonlinear)  BGK  type  model 
[87].  We  expect  these  results  to  have  an  impact  in  other  areas  of  applied  analysis. 

(ii) :  In  general  the  existence  theory  is  established  for  (a  system  of)  partial  (or  integro-partial) 
differential  equations  given  in  the  same  domain.  Our  models  (a:)  and  ((3)  consist  of  kinetic  equations 
and  their  hydrodynamic  limits  used  on  their  respective  domains  of  validity  and  coupled  through 
their  boundary  conditions.  We  thus  have  new  set  of  integro-differential  equations  with  new  types 
of  boundary  conditions.  As  far  as  we  know  no  such  types  of  problems  appeared  in  the  litterature 
before  (by  other  authors).  To  see  what  difficulties  we  must  face  to  handle  the  mathematical  theory 
of  these  novel  models,  assume  for  example,  that  we  apply  a  fixed  point  iterative  method  to  uncouple 
the  models.  Then  we  observe  that  we  must  deal  with  the  kinetic  equations  (resp.  the  hydrodynamic 
equations)  in  bounded  domains  with  new  types  of  boundary  conditions.  To  prove  the  existence  of  a 
fixed  point,  we  need' to  derive  appropriate  estimates  (particularly  at  the  interfaces  and  boundaries). 
It  is  not  straightforward  to  combine  the  kinetic  theory  with  the  hydrodynamic  theory  to  obtain  such 
estimates.  The  Sobolev  spaces  used  in  these  two  theories,  for  instance,  can  be  (and  are  indeed)  very 
different.  “Similar”  difficulties  arise  when  dealing  with  the  regularity,  positivity,  and  asymptotic 
behaviour  of  the  solutions  to  these  novel  models.  Therefore,  it  is  clear  that  new  methods  and  ideas 
have  to  be  introduced  and  developed  in  order  to  study  these  models.  Such  methods  have  been 
introduced  by  the  author  in  [77,  78,  82,  85,  80,  81].  Using  these  methods,  the  author  was  able  to 
establish  the  existence  theory  for  the  models  (a)  and  (f3)  for  the  case  of  Carleman  model,  linearized 
Boltzmann  model,  and  a  nonlinear  BGK  type  model  [80,  81,  87].  We  are  currently  working  on  the 
extension  of  these  results  to  more  general  models. 

(iii) :  The  asymptotic  analysis  and  trend  to  equilibrium  for  microscopic  models  is  in  general 
very  difficult  topic.  There  are  only  very  few  results  for  realistic  models.  For  the  models  (a)  and 
(/3),  the  situation  is  even  more  difficult  since  here  we  are  dealing  with  coupled  problems  of  two 
different  nature:  one  is  microscopic  and  the  other  is  macroscopic.  In  [77,  78,  82,  85,  80,  81]  we 
have  introduced  and  developed  a  general  method  that  led  to  the  asymptotic  analysis  and  trend 
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to  equilibrium  for  coupled  microscopic-microscopic  models  and  for  the  models  (o)  and  (/?)  for 
Carleman  model  and  a  linearized  Boltzmann  model.  This  method  also  allows  the  obtention  of  the 
rate  of  decay.  We  shall  use  these  ideas  to  obtain  further  results  on  the  asymptotic  and  trend  to 
equilibrium  models  (a)  and  (/?)  for  a  (nonlinear)  BGK  type  model. 

(iv):  This  part  is  even  more  difficult  since  to  obtain  any  analysis  of  the  approximations  of 
the  models  (a)  and  (f3)  by  the  transmission  time,  marching  algorithm,  we  need  to  know  detailed 
information  regarding  Parts  (ii)  and  (iii)  (the  continuous  problems).  Again  the  methods  introduced 
and  developed  by  the  author  [77,  78,  85,  82,  86,  80,  81]  were  shown  to  play  central  roles  in  the 
development  of  the  mathematical  foundations  of  the  approximations  of  the  models  (a)  and  (/?)  by 
the  transmission  time  marching  algorithm.  This  theory  is  obtained  for  the  case  of  Carleman  and 
linearized  Boltzmann  models  [78,  86,  80,  81].  Further  developments  are  underway. 


15 


References 

[1]  Andies,  P.  Le  Tallec,  P.,  Perlat,  J.  P.  and  Perthame,  B.,  The  Guassian  BGK  model  of  Boltzmann 
equation  with  small  Prandtl  numbers.  European  J.  Mechanics  (B  fluids),  813—830  (2000). 

[2]  Arkeryd,  L.  and  Cercignani,  C.,  A  global  existence  theorem  for  the  initial  boundary  value  problem 
for  the  Boltzmann  equation  whenthe  boundaries  are  not  isothermal,  Arch  Rat  Mech  Anal,  125, 
271-287,  (1993). 

[3]  Arlotti,  L.  and  Bellomo,  N.,  On  the  Cauchy  problem  for  the  nonlinear  Boltzmann  equation,  in: 
Lecture  Notes  on  the  Mathematical  Theory  of  the  Boltzmann  Equation,  World  Sci,  (1995). 

[4]  Babovsky,  H.,  On  a  simulation  schemefor  the  Boltzmann  equation,  Math  Meth  Appl  Sci,  8, 
223-233,  (1986) 

[5]  Babovsky,  H.,  A  convergence  proof  for  Nanbu’s  Boltzmann  scheme,  Europ  J  mech  Fluids,  8, 
41-55,  (1989). 

[6]  Babovsky,  H.  and  Illner  R.,  A  convergence  proof  for  Nanbu’s  simulation  method  for  the  full 
Boltzmann  equation,  SIAM  J  Num  Anal,  26,  45-65,  (1989). 

[7]  Bardos,  C.,  Golse,  F.,  and  Levermore,  D.,  Fluid  dynamics  limits  of  kinetic  equations  I:  Formal 
derivation  J.  Statistic  Phys,  63,  323-344,  (1991). 

[8]  Bardos,  C.,  Golse,  F.,  and  Levermore,  D.,  Fluid  dynamics  limits  of  kinetic  equations  7/Commun 
Pure  Appl  Math,  667-753,  (1993) 

[9]  Bellomo,  N.,  Palczewski,  A.  and  Toscani,  G.,  Mathematical  Topics  in  Nonlinear  Kinetic  Theory, 
World  Sci,  (1988). 

[10]  P.L.  Bhatnagar,  E.P.  Gross  and  M.  Krook,  Model  for  collision  processes  in  gases,  I.  Small 
amplitute  processes  in  charged  and  neutral  one-component  systems,  Phys.  Rev.  94,  511-524, 
(1954). 

[11]  Bourgat,  J.F.,  te  Tallec,  P.,  Mallinger,  F.,  Qiu,  Y.  and  Tidriri,  M.D.,  Numerical  coupling 
of  Boltzmann  and  Navier-Stokes,l.V.T.A.M.  (International  Union  of  Theoretical  and  Applied 
Mechanics)  Conference.  Marseille,  France  Sept.  1-4  1992. 

[12]  Bourgat,  J.F.,  Le  Tallec,  P.  and  Tidriri,  M.D.,  Coupling  Boltzmann  and  Navier-Stokes  Equa¬ 
tions  by  Friction  J.  Comp.  Phys.  127,  227-245  (1996). 

[13]  Brun,  R.,  Transport  et  relaxation  dans  les  ecoulements  gazeux ,  Masson,  Paris,  1986. 

[14]  Busoni,  G.  and  Palczewski  A,  Stationary  Boltzmann  equation  near  equilibrium,  M3AS:  Math 
Models  and  Methods  Appl  Sci,  3,  395-440,  (1993). 

[15]  Caflish,  R.  ,  The  fluid  dynamical  limit  of  the  nonlinear  Boltzmann  equation  Commun  Pure 
Appl  Math  33,  651-666  (1980) 

[16]  Cercignani,  C,  Theory  and  application  of  of  the  Boltzmann  equation,  Springer,  (1988). 

[17]  Chapman  S.,  The  kinetic  theory  of  simple  and  composite  gases:  viscosity,  thermal  conduction 
and  diffusion,  Proc.  Roy.  Soc.  (London)  A93  (1916/1917),  1—20. 


16 


[18]  Comeaux,  K.A.,  Chapman,  D.R.  and  MacCormack,  R.W.,  An  analysis  of  the  Burnett  equations 
based  on  second  law  of  Thermodynamics,  AIAA  Paper  95-0415,  January  1995. 

[19]  Coquel,  F.  and  Liou.  M.S.,  Hybrid  Upwind  Splitting  by  a  Field  by  Field  Decomposition,  NASA 
Lewis  Research  Center  Report. 

[20]  Chapman,  S.  and  Cowling,  T.G.,  The  mathematical  theory  of  non-uniform  gases  (4th  edn). 
Cambridge  University  Press  (1990). 

[21]  Darrozes  J  and  Guiraud  J,  Generalization  formelle  du  theoreme  H  en  presence  de  parois, 
Comp  Rend  Acad  Sci  Paris,  262A,  1368-1371,(1966). 

[22]  Desvillettes,  L.  and  Mischler,  S.,  Work  in  progress 

[23]  DiPerna,  R.  and  Lions,  P.L.  ,  On  the  Cauchy  problem  for  the  Boltzmann  -.Global  existence  and 
weak  stability  results,  Annals  of  Math,  130.  1189-1214,  (1990) 

[24]  DiPerna,  R.  and  Lions,  P.L.  ,  Global  solutions  of  the  Boltzmann  and  the  entropy  inequality, 
Arch  Rat  Mech  Anal,  114,  47-55,  (1991). 

[25]  DiPerna,  R.,  Measure-Valued  Solutions  to  Conservation  Laws,  Arch  Rat  Mech  Anal,  88,  223- 
270,  (1985). 

[26]  Enskog  D.,  Kinetische  Theorie  der  Vorgange  in  massig  verdunntent  Gasen,  1,  in  Allgemeiner 
Teil,  Almqvist  &  Wiksell,  Uppsala,  1917. 

[27]  Gilbarg,  D.  and  Trudinger,  N.S.,  Elliptic  partial  differential  equations  of  second  order ,  Springer 
Verlag  1977 

[28]  Golse,  F,  Lions,  P.L.,  Perthame,  B  and  Sentis,  R,  The  regularity  of  the  moments  of  the  solutions 
of  a  transport  equation,  J.  Funct  Anal.  76,  110-125,  (1988). 

[29]  Grad  H.,  Principles  of  the  kinetic  theory  of  gases,  In  Handbuch  der  Physik  (S.  Flugge  ed.), 
Vol.  XII,  Springer,  Berlin,  1958. 

[30]  Grad,  H.,  On  the  kinetic  theory  of  rarefied  gases,  Comm.  Pure  Appl.  Math.  2  331-407  (1947). 

[31]  Holway,  L.  H.,  Kinetic  theory  of  shock  structure  using  an  ellipsoidal  distribution  function. 
In  Rarefied  Gas  Dynamics,  Vol.  I  (Proc.  Fourth  Internat.  Sympos.,  Univ.  Toronto,  1964),  pp 
193-215.  Academic  Press,  New  York,  1966. 

[32]  Groth,  C.P.T.  and  Levermore,  C.D.  ,  Beyond  the  N avier-Stokes  approximation:  Transport 
corrections  to  the  Gaussian  closure,  in  preparation  (1996). 

[33]  Guiraud  J.P.  ,  Problemes  aux  limites  pour  I’equation  de  Boltzmann,  J.  Mecanique,  9,  443-490, 
(1970). 

[34]  Hamdache,  K.,  Initial-boundary  value  problems  for  the  Boltzmann  equation:  Global  existence 
of  weak  solutions  Arch  Rat  Mech  Anal,  bfll9309-353,  (1992). 

[35]  Hilbert  D .  Begrundung  der  kinetischen  Gastheorie.  Mathematische  Annalen  72  (1916/1917), 
562-577. 

[36]  M.  Kogan,  Rarefied  Gas  Dynamics ,  Plenum,  (1969). 


17 


[37]  Lions,  P.L.,  Compactness  in  Boltzmann  equation  via  Fourier  integrals  operator  and  applications 

I,  ,  J  Math  Kyoto  Univ,  bf34  391-427,  (1994). 

[38]  Lions,  P.L.,  Compactness  in  Boltzmann  equation  via  Fourier  integrals  operator  and  applications 

II,  ,  J  Math  Kyoto  Univ,  bf34  391-427,  (1994). 

[39]  Lions,  P.L.,  Compactness  in  Boltzmann  equation  via  Fourier  integrals  operator  and  applications 

III,  ,  J  Math  Kyoto  Univ,  bf34  391-427,  (1994). 

[40]  Lions,  P.L.,  Conditions  at  infinity  for  Boltzmann’s  equation,  Comraun  Partial  Diff  Equations, 
bfl9  335-367,  (1994). 

[41]  Lions,  P.L.,  Pertharae,  B.,  and  Tadmor,  E.,  A  kinetic  formulation  of  multidimentional  scalar 
conservation  laws  and  related  equations,  J .  Amer.  Math.  Soc.,  7  (1994)  169—191. 

[42]  Lions,  P.L.,  Perthame,  B.,  and  Tadmor,  E.,  Kinetic  formulation  of  the  isentropic  Gas  dynamics 
and  the  p-systems.  Commun.  Math.  Phys.  163  (1994) 

[43]  Lions,  P.L.,  Perthame,  B.,  and  Souganidis,  P.E.,  Existence  and  stability  of  entropy  solutions 
for  the  hyperbolic  systems  of  isentropic  gas  dynamics  in  Eulerian  and  Lagrangian  coordinates , 
(1996)  599-638. 

[44]  Mallet,  M.  A  finite  element  method  for  computational  fluid  dynamics,  Ph.D.  Thesis,  Stanford 
University,  (1985). 

[45]  Morse,  T.F.,  Kinetic  model  for  gases  with  internal  degrees  of  freedom,  Phys.  Fluids,  bf7  159-169 
(1964) 

[46]  Nanbu,  K.,  Direct  simulation  scheme  derived  from  the  Bolztmann  equation,  Japan  J  of  Phys, 
bf49  2042-2049,  (1980) 

[47]  Nanbu,  K.,  Stochastic  solution  method  of  the  Master  equation  and  the  model  Boltzmann  equa¬ 
tion,  J  Phys  Soc  Japan,  bf52  2654-2658,  (1983) 

[48]  Nanbu,  K.,  Stochastic  solution  method  of  Boltzmann  equation  I,  Memoirs  of  the  Institute  of 
Fluid  Sciences,  Tohoku  Univ,  bf3  47-93,  (1992). 

[49]  Legge  H,  Heat  transfer  and  forces  on  a  LIF  single- cristal  disc  in  hypersonic  flow,  DLR  Report 
No  222-92  A  14,(1992). 

[50]  Legge  H,  Force  and  heat  transfer  on  a  delta  wing  from  continuum  to  freem  molecular  flow,  in: 
Rarefied  Gas  Dynamics,  Am  Inst  Aeronaut  astronaut,  (to  appear). 

[51]  Le  Tallec,  P.  and  Tidriri,  M.D.,  Kinetic  upgrade  of  Navier-Stokes  equations  for  transitional 
flows.  Rapport  de  contrat,  6465/91  ATP  DPH.  October  1992. 

[52]  Le  Tallec,  P.  and  Tidriri,  M.D.,  Kinetic  upgrade  of  Navier-Stokes  equations  for  Transitional 
Flows,  Rapport  de  contrat  HERMES  R/Q  6465/91,  INRIA,  fevrier  1993. 

[53]  Le  Tallec,  P.  and  Tidriri,  M.D.,  Analysis  of  the  Time  Marching  Algorithm,  INRIA  Research 
Report  RR-2435,  Oct.  1994. 

[54]  Le  Tallec,  P.  and  Tidriri,  M.D.,  Convergence  Analysis  of  Domain  Decomposition  algorithms 
with  full  overlapping  for  the  advection-diffusion  problems,  ICASE-TR  No.  96-37,  June,  1996. 


18 


[55]  Le  Tallec,  P.  and  Tidriri,  M.D.,  Analysis  of  the  explicit  Time  Marching  Algorithm,  ICASE 
Report  No.  96-45,  August,  1996. 

[56]  Le  Tallec,  P.  and  Tidriri,  M.D.,  Convergence  of  domain  decomposition  algorithms  with  full 
overlapping  for  the  advection-diffusion  problems.  Accepted  for  publication  in  Math.  Comp. 

[57]  Le  Tallec,  P.  and  Perlat,  J.-Ph.,  Analyse  numerique  des  systemes  aux  moments  de  Levermore , 
Rapport  de  recherche  INRIA  Octobre  1996. 

[58]  Levermore,  C.D.  ,  Moment  Closure  Hierarchies  for  Kinetic  Theories,  J.  Stat.  Phys.  83  1021- 
1065, (1996). 

[59]  Murat,  F.  ,  Compacite  par  compensation,  Ann.  Scuola.  Norm.  Sup.  Pisa  Sci.  Fis.  Math.  bf5 
489-507,  (1978). 

[60]  Nishida,  T.  ,  Fluid  dynamical  limit  of  the  nonlinear  Boltzmann  equation  to  the  level  of  the 
incompressible  Euler  equation,  Commun.  Math.  Phy.,  bf61  119-148  (1978). 

[61]  Palczewski,  A.,  Spectral  Properties  of  the  space  inhomogeneous  linearized  Boltzmann  operator, 
Transport  Theory  Statist  Phys  bfl3  409-430  (1984) 

[62]  Perthame,  B.,  Global  existence  to  the  BGK  model  of  the  Boltzmann  equation,  J.  Differential 
Equations,  bf82  191-205  (1989). 

[63]  Perthame,  B.,  Kinetic  equations  and  hyperbolic  systems  of  conservation  laws,  in  proceedings 
of  the  international  congress  of  mathematicians,  Zurich,  Switzerland  1994. 

[64]  Perthame,  B.,  and  Pulvirenti,  M.,  On  some  large  systems  of  random  particles  which  approxi¬ 
mate  scalar  conservation  laws,  Asymptotic  Anal.  10  (3)  (1995),  263-278. 

[65]  Perthame,  B.,  and  Pulvirenti,  M.,  Weighted  L°°  bounds  and  uniqueness  for  the  Boltzmann 
BGK  model,  Arch.  Rational  Mech.  Anal.  125  (1993)  289-295. 

[66]  Pulvirenti,  M.,  Wagner,  W.and  Zavelani  Rossi,  M.,  Convergence  of  particle  schemes  for  the 
Boltzmann  equation,  Eur  J  Mech  B/Fluids,  bfl3  339-351,  (1994). 

[67]  Ringeissen,  E,  Thesis,  University  of  Paris  VII,  1991. 

[68]  Rostand,  Ph.  and  Stoufflet,  B.,  Finite  volume  Galerkin  methods  for  viscous  gas  dynamics,  Rap¬ 
port  de  recherche  INRIA  no  863,  Juillet  1988. 

[69]  Shakib,  F.,  Finite  element  analysis  of  the  compressible  Euler  and  Navier-Stokes  equations 
Ph.D.  Thesis,  Stanford  University,  (1988). 

[70]  Shizuta,  Y.,  On  the  classical  solutions  of  the  Boltzmann  equation,  Commun  Pure  Appl  Math 
bf36  705-754  (1983) 

[71]  Sideris,  T.,  Formation  of  singularities  in  three  dimensional  compressible  fluids,  Commun  Math. 
Phy.  bflOl  475-485  (1985). 

[72]  Tartar,  L.,  Une  nouvelle  methode  de  resolution  d’equations  aux  derivees  partielles  nonlineaires, 
in  Lecture  Notes  in  Mathematics,  bf665  Springer,  228-241  (1977). 

[73]  Tidriri,  M.D.,  Couplage  d’ approximations  et  de  modeles  de  types  differents  dans  le  calcul 
d’ecoulements  externes,  PhD  thesis,  Universite  de  Paris  Dauphine,  1992. 


19 


[74]  Tidriri,  M.D.,  Domain  Decomposition  for  Incompatible  nonlinear  models,  INRIA  Research 
Report  RR-2378,  October  1994. 

[75]  Tidriri,  M.D.,  Domain  Decompositions  for  Compressible  Navier-Stokes  Equations  with  different 
discretizations  and  formulations,  J.  Comp.  Phys.  119,  271-282  (1995). 

[76]  Tidriri,  M.D.,  Mathematical  Analysis  of  the  Navier-Stokes  Equations  with  Non  Standard 
Boundary  Conditions,  ICASE-TR  No.  95-54,  July  1995. 

[77]  Tidriri,  M.D.,  Asymptotic  analysis  of  a  coupled  system  of  kinetic  equations.  C.  R.  Acad.  Sci. 
Paris,  t.  328,  Serie  I  Math.,  637-642,  1999. 

[78]  Tidriri,  M.D.,  Numerical  analysis  of  coupling  for  a  kinetic  equation.  Mathematical  Modelling 
and  Numerical  Analysis  M2AN,  Vol.  33,  No  6,  1999,  p.  1121-1134. 

[79]  Tidriri,  M.D.,  Hybrid  finite  element/finite  volume  approximation  to  conservation  laws:  Exis¬ 
tence  theory.  Comm.  Appl.  Anal.,  Vol.  5,  2001,  no.  1,  pp.  91-101. 

[80]  Tidriri,  M.D.,  Rigorous  derivation  and  analysis  of  coupling  of  kinetic  equations  and  their  hy¬ 
drodynamics  limits  for  a  simplified  Boltzmann  model.  Journal  of  Statistical  Physics,  Vol.  10Jt. 
Nos.  1/2  (2001),  pp.  257-292. 

[81]  Tidriri,  M.D.,  New  models  for  the  solution  of  intermediate  regimes  in  transport  theory  and 
radiative  transfer:  existence  theory,  positivity,  asymptotic  analysis,  and  approximations.  Journal 
of  Statistical  Physics,  Vol.  1  Of.  Nos.  1/2  (2001),  pp.  293-328. 

[82]  Tidriri  M.  D.,  A  novel  class  of  multiscale  models  in  mathematical  physics.  Journal  of  nonlinear 
analysis,  in  press,  2001. 

[83]  Tidriri  M.  D.,  An  alternative  model  to  boundary  and/or  interface  layers  correctors  in  semi¬ 
conductor  theory:  derivation  and  existence  theory  C.  R.  Acad.  Sci.  Paris,  t.  332,  No.  3.,  Serie  I 
Math.,  pp.  271-276. 

[84]  Tidriri  M.  D.,  Alternative  models  to  boundary  and/or  interface  layers  correctors  transport 
theory  and  radiative  transfer.  C.  R.  Acad.  Sci.  Paris,  Serie  I  (Math.)  333  (3)  199-204  (2001) 

[85]  Tidriri  M.  D.,'  Asymptotic  analysis  of  the  coupling  of  transport  equations  and  their  diffusion 
approximations.  C.  R.  Acad.  Sci.  Paris,  t.  330,  Serie  I  (Math.),  p.  1073-1078,  2000. 

[86]  Tidriri  M.  D.,  Applications  of  the  transmission  time  marching  algorithm  to  the  coupling  of 
transport  equations  and  their  diffusion  approximations.  C.  R.  Acad.  Sci.  Paris,  Serie  I  (Math.) 
333  (4)  (2001). 

[87]  Tidriri,  M.D.,  Analysis  of  a  coupled  system  of  kinetic  equations  and  conservation  laws.  Ac¬ 
cepted  in  C.  R.  Acad.  Sci.  Paris,  Serie  I  (Math.) 

[88]  Tidriri,  M.D.,  Work  in  progress. 

[89]  Zhong,  X.,  MacCormack,  R.W.  and  Chapman,  D.R.,  AIAA  J.  bf31  6  (1993) 


20 


\ 


02  0.03  0.04  0.05  0.06  0.07  0.08 

x-coordinates 


ure  1:  Skin  Friction  coefficient  for  Kn 

10“3 


0.02  0.03 


0.04  0.05  0.06 

x-coondinates 


0.07  0.08 


Figure  3:  Pressure  coefficient  for  Kn  =  8.310 
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Figure  2:  Heat  Flux  coefficient  for  Kn  —  8.310  3 
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Figure  4:  Density  Cross  Section  for  Kn 
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